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CONTRACTOR REPORT 


MISSILE AERODYNAMICS FOR ASCENT AND RE-ENTRY 

INTRODUCTION 


This document presents equations for aerodynamic forces and moments to be used in 6-DOF 
missile simulations such as MAVERIC. The Missile coordinate frame (M frame) and a frame 
parallel to the M frame were used for aerodynamics to allow convenient simulation of either the 
ascent phase of flight or a tumbling re-entry. The missile configuration chosen as an example is a 
cylinder with fixed fins and a nose cone. The equations include both the static aerodynamic 
coefficients and the aerodynamic damping derivatives. The inclusion of aerodynamic damping is 
essential for simulating a tumbling re-entry in which large angular rates may occur. The missile 
aerodynamic reference area and length are the cylinder cross-sectional area (S re j) and diameter 
(D re f), respectively. 

Appendix A presents a mathematical equivalency for damping derivatives to give insight into 
aerodynamic damping. Appendix B provides additional insight into aerodynamic damping by 
presenting an analytical technique for calculating the damping moment of a cylinder at its center- 
of-mass in cross-flow conditions. 



STATIC COEFFICIENT MODEE 


The static coefficients in this document are dimensionless and are defined in the M frame. The M 
frame, shown in Figure 1 , can greatly simplify the missile aerodynamics, as will be discussed 
later in the document. The M frame naturally fits the situation of large angular excursions 
experienced by a tumbling missile and allows an analyst to more easily visualize the missile 
rotational dynamics. The M-frame origin is at the moment reference point (MRP) on the missile 
centerline with the X-axis pointing forward along the centerline. For purposes of illustrating the 
M frame, the wind-relative velocity vector, V R , has been translated from the missile center-of- 
mass to the MRP. The Y-axis of the M frame is pointed in the direction of V R x l , where i is 
a unit vector along the X-axis of the M frame. In other words, V R (after being translated to the 
MRP) remains in the X-Z plane of the M frame regardless of the roll position of the missile. The 
total angle-of-attack, cc tot , shown in Figure 1, is measured in the X-Z plane of the M frame and 
has a range of 0° to +180°. 




Figure 1. M frame, P frame, and static coefficients. 
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The six static coefficients, shown in Figure 1, are listed below. The final subscript “m” denotes 
the M frame. 


Q Am 

axial force coefficient 

Cym 

side force coefficient 

(-Nm “ 

normal force coefficient 

Clm — 

rolling moment coefficient 

r — 

° mm 

pitching moment coefficient 

r — 

'-‘yawm 

yawing moment coefficient 


The axial force coefficient, C Am , and the nonnal force coefficient, C Nrn , are positive in a 
direction opposite to that of the X-axis and Z-axis, respectively. The three moment coefficients 
are positive by the right-hand rule. 

C Am , C Nrn , and C mm are called the pitch plane (or longitudinal) coefficients. C Ym , C lm , and 
Cyawm are called the lateral-directional (or lateral) coefficients and can be set equal to zero in the 
M frame for a missile that is axisymmetric about the centerline. The reason that the lateral- 
directional coefficients can be set equal to zero for an axisymmetric missile (leaving only the 
axial force, nonnal force, and pitching moment) is that the X-Z plane of the M frame is a plane 
of symmetry containing the V R vector. 

Figure 1 also shows a body-fixed P frame, which is called the MRP frame, with origin at the 
MRP and with axes parallel to those of a traditional Body frame (B frame, not shown). The 
origin of the B frame is at the missile center-of-mass with the X-axis forward, the Y-axis 
starboard, and the Z-axis completing the right-handed system. 

The aerodynamic roll angle, <p A , is the rotation angle between the M and P frames and will be 
defined mathematically in the next paragraph. The P-frame angle-of-attack, a p , and sideslip 
angle, /? , in Figure 1, are not used in the formulation of static coefficients in the M frame. 
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Static coefficients for a missile are, in general, a function of a tot , cp A , and Mach number. For a 
missile that is axisymmetric about the centerline, the static coefficients in the M frame are only a 
function of a tot and Mach, because the aerodynamic properties do not change with cp A . The 
equations for these three variables are: 



^ = t ™- 1 (^) 

Vr 

MACH = 

* sound. 

where V Rx , V Ry , V Rz are components of V R in the B frame, translated to the MRP 

V sound. is the speed of sound at the missile’s altitude 

a tot has a range of 0° to +180° and is never negative 

(p A has a range of - 1 80° to + 1 80° 


( 2 ) 

( 3 ) 


It is imperative to note that singularities exist in the <p A equation at a tot = 0° and 180°, and these 
must be handled by specifying a value of (p A (such as 0°) for these two values of a tot . 


If aerodynamic uncertainties are to be modeled in a simulation, the dispersions should be applied 
to the static coefficients in the M frame. The dispersed static coefficients will then be used to 
calculate aerodynamic forces and moments on the missile in the M frame as described later in 
this document. 
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AERODYNAMIC DAMPING MODEL 


The missile aerodynamic damping model is based on traditional damping derivatives that are 
dimensionless and are used to calculate aerodynamic moments caused by angular rates of the 
missile. See Appendix A for mathematical insight into damping derivatives. Missile damping 
derivatives are typically negative and, therefore, tend to reduce the angular rates. 

Damping derivatives should be developed at (or near) the missile center-of-mass rather than at 
the MRP. See Appendix B for an explanation of the effects of longitudinal center-of-mass 
location upon damping. Damping derivatives should also be developed in a coordinate frame 
parallel to the M frame to properly separate the M-frame pitch and yaw damping effects. The 
coordinate frame that will be used for calculating missile damping derivatives is parallel to the 
M frame and has its origin on the centerline at the same X-axis location as the center-of-mass. In 
using such a coordinate frame on the centerline, the effects upon damping of small lateral center- 
of-mass offsets from the centerline are ignored. 

Damping derivatives represent pure couples at the center-of-mass, and these couples may be 
applied at any other point. After being calculated on the centerline at the same X-axis location as 
the center-of-mass, the damping derivatives will be applied at the MRP in the M frame. Damping 
derivatives will be given names as though they were originally established in the M frame. 
Damping derivatives in this document are Ci pm , C mqm , and C yawrm for roll, pitch, and yaw, 
respectively. These damping derivatives are positive in the M frame by the right-hand rule, and 
their mathematical definitions are: 


^Ipm 


dC, 




Pm A 


ref 


2V, 


R 


c 


mqm 


dC 


mm 


q m D 


ref 


2V , 


R 


c. 


yawrm 


dC 


yawm 


r m D r ef 


2V, 


R 


( 4 ) 

( 5 ) 

( 6 ) 


The dimensionless tenns in the denominators of equations (4), (5), and (6) do seem a bit strange, 
but their mathematical basis can be found in a close examination of the damping equations in 
Appendix A. 
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These damping derivatives use M-frame angular rates (p m , q m , r m ), which are obtained by 
transforming the body rates from the B frame to the M frame by a negative rotation ( -<p A ) about 
the centerline. For more complicated vehicle configurations such as airplanes, aerodynamicists 
sometimes use additional damping moment terms that are not discussed in this document. There 
also might be a change in an aerodynamic force due to an angular rate, such as the change in 
nonnal force caused by a pitch angular rate, but this type force will not be considered at this 
time. 


At (Xtot = 0°, the pitch damping derivative for the missile is assumed to be equal to the yaw 
damping derivative (also applies for a, ot =180°) because the missile shape approximates a 
vehicle that is axisymmetric about the longitudinal axis. The yaw damping derivative for the 

missile at a tot = 90° (cross-flow conditions) is assumed to be much less than the corresponding 
pitch damping derivative because of the difference in the orientations of the pitch and yaw axes 
with respect to the wind-relative velocity vector. 


If aerodynamic uncertainties are to be modeled in a simulation, the dispersions should be applied 
to the damping derivatives in the M frame. The dispersed damping derivatives should then be 


converted into the form of dimensionless damping moment coefficients ( Ci mcL , C mmd , C yawmd ) 

in the M frame by assuming a linear relationship between coefficients and angular rates: 


C Imd 

Pm D r ef 
2 V r 

Clpm 

( 7 ) 

Cmmd 

q m D re f 
2 V r 

■ r 

^ mqm 

( 8 ) 

f ' _ 

Cyawmd 

T m D re f 
2 V r 

r 

° yawrm 

( 9 ) 


The damping moment coefficients in equations (7), (8), and (9) can be interpreted as the changes 
in the rolling, pitching, and yawing moment coefficients caused by damping. Damping moment 
coefficients will be used with the six static coefficients to calculate the aerodynamic forces and 
moments. 
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AERODYNAMIC FORCES AND MOMENTS 


Aerodynamic forces and moments will first be calculated in the M frame at the MRP using the 
six static coefficients and three damping moment coefficients described in this document. The 
vector equations for the aerodynamic forces and moments in the M frame are: 


F = 

1 m 




M m = 


M X m\ I 

( Ci 

° lm 

+ 

C Imd 

M Ym ) — Q Sref D ref 

c 

^ mm 

+ 

Cmmd 

M Z J \ 

c 

y ° yawm 

+ 

Cyawmd. 


( 10 ) 


( 11 ) 


where Q is the missile dynamic pressure. These aerodynamic forces and moments must be 
transformed from the M frame to the P frame by a rotation about the centerline through the 
aerodynamic roll angle, (p A . These transfonnations, in matrix notation, are: 


Xp 

U Yp j = 

v 


0 

cos{(p A ) 
- sin((p A ) 


0 

sin(cp A ) 
cos(<p A ) 



'M xp 

Mv-n I = 


M 


l Yp 

Zp, 


0 

cos{(p A ) 


-sinifpA) cos((p A ) 



( 12 ) 


( 13 ) 


These P-frame aerodynamic forces and moments must finally be transfonned to the B frame at 
the missile center-of-mass, taking into account the location of the moment reference point 
(. MRP b ) relative to the center-of-mass. The vector equations for the aerodynamic forces and 
moments in the B frame are: 
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(M x i A / M xA /MRP^n /fx P \ 

M b — ( ^y& I — ( ] + I MRP Y b ) x J F Yp ) (15) 

\M zb J \M Zp ) \MRP zb ) \F Zp ) 

Carrying out the cross product in equation (15), the moment equations in the B frame become: 

M xb — M X p + MRP Yb F Z p — MRP zb F Y p ( 16 ) 

M Yb — M Y p + MRP zb F X p — MRP xb F Zp ( 17 ) 

M zb — M Z p + MRP xb F Y p — MRP Yb F Xp ( 18 ) 

The force and moment equations in this section should be used to incorporate the aerodynamic 
models in this document into a 6-DOF simulation. 
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APPENDIX A— 


MATHEMATICAL EQUIVALENCY FOR DAMPING DERIVATIVES 


Nomenclature 

Cf in normal force coefficient of fin based on fin projected area, Sfin 

C mqm pitch damping derivative 

— > 

Ffin normal force on fin 

M d damping moment caused by fin 

Mf in total pitching moment caused by fin 

Q d dynamic pressure for fin damping calculation 

Qf in dynamic pressure at fin 

q m pitch angular rate (radians/s) 

Vr magnitude of the wind-relative velocity vector at center-of-mass 

This Appendix provides insight into aerodynamic damping by showing that a pitch damping 
derivative, C mqrn , for a missile is equivalent to the negative of a nonnal force coefficient for a fin 
on the missile. The example used to develop the mathematics for this coefficient equivalency is a 
simple missile (not necessarily realistic) with a horizontal fin at the rear to provide all the pitch 
damping (Figure A.l). The distance from the fin to the missile center-of-mass will be set equal to 
the missile reference length, D re f, solely for the purpose of obtaining a final result without 
inserting additional constants. For this derivation the M, P, and B frames (not shown) are 
assumed to be coincident. All three frames have their origins on the centerline at the center-of- 
mass with their Z-axes pointed downward. The equations are simplified by assuming cross-flow 
conditions. 


D ref 



Figure A. 1 Example missile for damping derivative equivalency. 
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By inspection of Figure A.l, the dynamic pressure at the fin is: 


1 

Qfin ~ 2 P R + D r ef Qm + ^ref Qm) (A‘ 1 ) 

The magnitude of the third term in parentheses in equation (A.l) is typically less than 1% of the 
magnitude of the second term, so the third term will be ignored. The fin dynamic pressure now 
becomes: 


Qfin = 2 P ( V R + 2 V R D ref ) (A- 2 ) 

The fin nonnal force is simply: 

Ffin — Qfin C fin ^ fin (A- 2 ) 

where Sj in is the fin projected area and Cf in is the normal force coefficient of the fin based on 
the fin projected area. The total pitching moment caused by the fin is the negative of the fin 
nonnal force times the moment arm, D re f, as follows: 

Mfin — — Qfin Cfin $ fin D ref A- 4 ) 

It will now be assumed that the fin projected area, Sf in , is 1/4 of the missile reference area, S re f . 
(This assumption permits the final result to be obtained without inserting additional constants). 
Plugging this 1/4 fin area relationship into equation (A.4) produces: 



1 

^ Qfin Cfin $ ref D re f 


(A 5) 


Going back to equation (A.l) for the dynamic pressure at the fin, the first tenn inside the 
parentheses is for the static moment, so only the second tenn is needed for damping calculations. 
The dynamic pressure to be used for fin damping is therefore: 

Qd ~ P Vr Dref Qm (A- 


Substituting the right side of equation (A. 6) for the Qf in tenn in equation (A. 5) produces the fin 
damping moment as follows: 


M d = 


1 

^ P Vr Qm Cfin ^ref 



(A 7) 


10 



The fin damping moment will now be expressed in a generic form using the following traditional 
damping moment coefficient equation: 


M d = 


2 

— P Vr S re f D r ef 


C- 


mmd 


(A 8) 


where C mrnd is the pitch damping moment coefficient, which represents the change in the 
pitching moment coefficient caused by fin damping in this example. Combining equations (A. 7) 
and (A. 8) to eliminate M d produces the following equation that shows a linear relationship 
between C rnmd and q m for the simple example chosen: 


c 


Pm^ref 


mmd 


2V f 


C 


R 


fin 


(A 9) 


C-mmd can also be expressed, in a general sense, as the following partial differential equation: 

r _ d C m m 

mmd 3 “ ' Pm 

OC[m 

Combining equations (A. 9) and (A. 10) to eliminate C mmd , and then solving for the normal force 
coefficient of the fin, Cf in , produces the following relationship: 


(A. 10) 



dC 


mm 


Pm A 


ref 


2 V, 


R 


(All) 


The traditional equation that defines the pitch damping derivative, C mqm , is: 


c 


dC 


mm 


mqm 


Pm D 


ref 


2 V, 


R 


(A. 12) 


The identical expressions on the right side of equations (A. 11) and (A. 12) illustrate that for 
purposes of missile damping, a pitch damping derivative is equivalent to the negative of a nonnal 
force coefficient for a horizontal fin. 


11 



APPENDIX B— 


PITCH DAMPING MOMENT FOR A CROSS-FLOW CYLINDER 


Nomenclature 

C cylinder cross-flow drag coefficient (The reference area for C is always the 

cylinder cross-flow projected area being analyzed, whether that area is an 
incremental area or the total area; C is assumed to not vary with the size of 
the cylinder or the size of any incremental projected area) 
eg distance from center-of-mass to forward end of cylinder, expressed as a fraction 

of cylinder length (forward end is right side of cylinder in figure B.l) 

D cylinder diameter (total projected area in cross-flow conditions = L D) 

dx incremental distance along cylinder length (incremental projected area = dx D) 

Fj incremental aerodynamic force on incremental projected area (Fj = Q C dx D) 

L cylinder length 

M d total damping moment about center-of-mass 

Mi incremental moment about center-of-mass 

Q dynamic pressure at incremental projected area 

q m pitch angular rate (radians/s) 

V R magnitude of the wind-relative velocity vector at center-of-mass 

X distance from center-of-mass to incremental projected area (positive forward) 

This Appendix provides additional insight into aerodynamic damping by presenting a derivation 
of the pitch damping moment for a cylinder in cross-flow conditions (Figure B.l). The effects of 
the longitudinal location of the center-of-mass are included and are very important. The cylinder 
is assumed to have the aerodynamic properties of an infinite cylinder (no end flow). For this 
derivation, the M, P and B frames (not shown) are assumed to be coincident. All three frames 
have their origins on the centerline at the center-of-mass with their Z-axes pointed downward. 


Mi,q r 


F ; 


d T 


Center of Mass 


(cg-L-L) 


dx 


(cg-L ) 


V. 


Figure B.l Cross-flow cylinder. 
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By inspection of Figure B.l, the dynamic pressure at the incremental projected area is: 

Q = 2 p( ~ v * - 2V ^ X + 9m* 2 ) (B - 1} 

Typically, the magnitude of the third term in parentheses in equation (B.l) is approximately 
1% of the magnitude of the second term, so the third term will be ignored. The dynamic 
pressure now becomes: 

<3 = \p(Vr ~ 2V R q m X) (B.2) 

The incremental aerodynamic force on the incremental projected area is: 

Ft = \p(Vi - 2 V R q m X) CDdx (B. 3) 

The incremental moment about the center-of-mass caused by the incremental aerodynamic 
force is: 

M,= ^p(Vl - 2V R q m X) C X D dx (B.4) 

The first tenn in parentheses of equation (B.4) is for calculating the incremental static moment 
and is not needed in this derivation. The second term, —2 V R q m X , is the term of interest in this 
example because it determines the incremental damping moment, M di , as follows: 

M di = -pV R q m C X 2 D dx (5. 5) 

Integrating equation (B.5) from the aft end to the forward end of the cylinder (from left to 
right in Figure B.l) produces the following equation for the total damping moment about the 
center-of-mass: 

M d = -\pV R qm C D L 3 (l - 3 eg + 3 eg 2 ) (5.6) 

Equation (B.6) for the total damping moment is the final result for this example, and it shows 
the significant sensitivity of pitch damping to the longitudinal center-of-mass location and to 
the cylinder length. 


13 



REPORT DOCUMENTATION PAGE 

Form Approved 
OMB No. 0704-0188 

The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and maintaining 
the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing 
this burden, to Department of Defense, Washington Headquarters Services, Directorate for Information Operation and Reports (0704-0188), 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302. 
Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to any penalty for failing to comply with a collection of information if it does not display a currently valid 
OMB control number. 

PLEASE DO NOT RETURN YOUR FORM TO THE ABOVE ADDRESS. 

1. REPORT DATE (DD-MM-YYYY) 
01-11-2012 

2. REPORT TYPE 

Contractor Report 

3. DATES COVERED (From - To) 

October 201 1-September 2012 

4. TITLE AND SUBTITLE 


5a. CONTRACT NUMBER 

NNM12AA41C 

Missile Aerodynamics for Ascent and Re-entry 

5b. GRANT NUMBER 



5c. PROGRAM ELEMENT NUMBER 

6. AUTHOR(S) 


5d. PROJECT NUMBER 

G.L. Watts* and J.W. McCarter** 


5e. TASK NUMBER 



5f. WORK UNIT NUMBER 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

George C. Marshall Space Flight Center, Huntsville, AL 35812 
Jacobs ESSSA Group, Huntsville, AL 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

M-1351 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 

10. SPONSORING/MONITOR’S ACRONYM(S) 

Washington, DC 20546-0001 


11. SPONSORING/MONITORING REPORT NUMBER 

NASA/CR — 20 1 2-2 1 7475 

12. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified-Unlimited 
Subject Category 02 

Availability: NASA CASI (443-757-5802) 

13. SUPPLEMENTARY NOTES 

Prepared for the Spacecraft & Vehicle Systems Dept., Engineering Directorate 

COTR: Kimball Y. Ibrahim 

*Oualis Corporation **DESE Research 

14. ABSTRACT 

Aerodynamic force and moment equations are developed for 6-DOF missile simulations of both the 
ascent phase of flight and a tumbling re-entry. The missile coordinate frame (M frame) and a frame paral- 
lel to the M frame were used for formulating the aerodynamic equations. The missile configuration chosen 


as an example is a cylinder with fixed fins and a nose cone. The equations include both the static aero- 
dynamic coefficients and the aerodynamic damping derivatives. The inclusion of aerodynamic damping 
is essential for simulating a tumbling re-entry. Appended information provides insight into aerodynamic 
damping. 


15. SUBJECT TERMS 


missile aerodynamics, aerodynamic damping, damping derivatives, missile coordinate frame, missile 
re-entry, 6 degrees-of-freedom 


16. SECURITY CLASSIFICATION OF: 

17. LIMITATION OF ABSTRACT 

18. NUMBER OF 

19a. NAME OF RESPONSIBLE PERSON 

a. REPORT 

b. ABSTRACT 

c. THIS PAGE 

UU 

PAGES 

STI Help Desk at email: help@sti.nasa.gov 

u 

u 

u 

28 

19b. TELEPHONE NUMBER (Include area code) 

STI Help Desk at: 443-757-5802 


Standard Form 298 (Rev. 8-98) 
Prescribed by ANSI Std. Z39-18 



National Aeronautics and 
Space Administration 
IS 20 

George C. Marshall Space Flight Center 

Huntsville, Alabama 35812 


